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Abstract 

It is shown that eightfold degeneracy in neutrino oscillations is easily seen by plotting constant 
probabilities in the (sin 2 28\^, l/s 2 ^) plane. Using this plot, we discuss how an additional long 
baseline measurement resolves degeneracies after the JPARC experiment measures the oscillation 
probabilities Piy^ — ► f e ) and Piy^ — ► v e ) at \Arri^i\L/4E = tt/2. By measuring -P(z^ — > v e ) 
or P{Vfj, — ► u e ), the sgn(Am3 1 ) ambiguity is resolved better at longer baselines and the 5 «-> 
it — 5 ambiguity is resolved better when |Am| 1 |L/4£' — vr/2| is larger. The #23 ambiguity may 
be resolved as a byproduct if — n\ is small and the CP phase 6 turns out to satisfy 

|cos(e) + I Am'^ 1 \L/4:E)\ ~ 1. It is pointed out that the low energy option (E ~lGeV) at the off- 
axis NuMI experiment may be useful in resolving these ambiguities. The v e — > v T channel offers a 
promising possibility which may potentially resolve all the ambiguities. 
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I. INTRODUCTION 



From the recent experiments on atmospheric [l| and solar |2j, and reactor [3|,y| neutrinos, 
we now know approximately the values of the mixing angles and the mass squared differences 
of the atmospheric and solar neutrino oscillations: 



(sin 2 20i2, Am^) ~ (0.8,7x10 5 eV 2 ) for the solar neutrino 
(sin 2 26*23, |Am 31 |) ~ (1.0,2 x 10~ 3 eV 2 ) for the atmospheric neutrino, 

where we use the standard parametrization of the MNS mixing matrix 

— i5 

C12C13 S12C13 s^e 

U = I -S12C23 - Ci 2 S 2 3Sl3e i<5 C12C23 - Si 2 S 23 Si 3 e lS S23C13 
■S12-S23 — Ci2C23Si3e* — C12S23 — Si2C23Si3e J C23C13 

and the case of Am^ > (Am^ < 0) corresponds to the normal (inverted) mass hierarchy, 
as is shown in Fig.0 In the three flavor framework of neutrino oscillations, the oscillation 
parameters which are still unknown to date are the third mixing angle #13, the sign of the 
mass squared difference Am^ of the atmospheric neutrino oscillation, and the CP phase 
5. It is expected that long baseline experiments in the future will determine these three 
quantities. 

Since the work of [6] , it has been known that even if the values of the oscillation probabil- 
ities Piy^ — > u e ) and Pip^ — > u e ) are exactly given we cannot determine uniquely the values 
of the oscillation parameters due to parameter degeneracies. There are three kinds of param- 
eter degeneracies: the intrinsic (#13, 8) degeneracy the degeneracy of AmJi <-> — Am|i 
0|, and the degeneracy of #23 «-> 7r/2 — #23 [H ©]. The intrinsic degeneracy is exact when 
Amli/ArnJi is exactly zero. The sgn(Am| 1 ) degeneracy is exact when AL is exactly zero, 
where A = \/2GpN e and L stand for the matter effect and the baseline, respectively (Gp is 
the Fermi constant and N e is the electron density in matter). The #23 degeneracy is exact 
when cos 2^23 is exactly zero. Each degeneracy gives a twofold solution, so in total we have 
an eightfold solution if all the degeneracies are exact. In this case prediction for physics is the 
same for all the degenerated solutions and there is no problem. However, these degeneracies 
are lifted slightly in long baseline experiments 1 , and there are in general eight different so- 
lutions 0. When we try to determine the oscillation parameters, ambiguities arise because 
the values of the oscillation parameters are slightly different for each solution. In particu- 
lar, this causes a serious problem in measurement of CP violation, which is expected to be 
small effect in the long baseline experiments, and we could mistake a fake effect due to the 
ambiguities for nonvanishing CP violation if we do not treat the ambiguities carefully. 

In the references 0, 0,111 in the past, various diagrams have been given to visualize how 
degeneracies are lifted in the parameter space. To see how the eightfold degeneracy is lifted, 
it is necessary for the plot to give eight different points for different eight solutions. An 



effort was made in 11] to visualize the eight different points by plotting the trajectories of 
constant probabilities in the (sin 2 2#i 3 , S23) plane. In the present paper we propose a plot 



1 cos 2#23 may be exactly zero, but the present atmospheric neutrino data |lo| allow the possibility of 
cos 26>23 7^ 0, so we will assume cos 2#23 7^ in general in the following discussions. 
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in the (sin 2 2#i 3 , l/s| 3 ) plane, which offers the simplest way to visualize how the eightfold 
degeneracy is lifted. As a byproduct, we show how the third measurement of — > v e , 
Vp, — > v e or v e — > u T resolves the ambiguities, after the JPARC experiment l3| measures 



the oscillation probabilities Piy^ — > u e ) and Piy^ — > u e ) at the oscillation maximum, i.e., at 
\Am 2 31 \L/4E = n/2. 

In the following discussions we assume that (Am 2 ,!, Am^ and 9i 2 are sufficiently precisely 
known. This is justified because the correlation between these parameters and the CP phase 
5 is not so strong in the case of JPARC [13], and we can safely ignore the uncertainty of 
these parameters to discuss the ambiguities in 6 due to parameter degeneracies. 



II. PLOTS IN THE (sin 2 26 13 , l/s| 3 ) PLANE 

As in Ref. let us discuss the ambiguities due to degeneracies step by step in the 

order {9 23 - tt/4 = 0, Am 2 21 = 0, A = 0) -> {9 23 - tt/4 ^ 0, Am 2 21 = 0, A = 0) -> (# 23 - vr/4 ^ 
0, Am 2 , ^ 0, A = 0) -> (9 23 - tt/4 ^ 0, Am 2 21 ^ 0, A ^ 0). 



A. cos 2<9 23 = 0, Am 2 2l /Am 2 3l = 0, AL = 

In this case the oscillation probabilities Piv^ — > f e ) and -P(^ — > v e ) are equal and are 
given by 

Piy^ -> u e ) = P(u lt -> v e ) = S23 sin 2 29 13 sin 2 A, 
where we have introduced the notation 

|Am 31 |L 



A 



AE 



To plot the line Piy^ — > u e ) = Piy^ — > z/ e )=const. in the (sin 2 2#i 3 , l/s| 3 ) plane, let us 
introduce the variables 

X = sin 2 2# 13 , 

r = ±- 

s 23 



Then 



give a straight line 



P = s 2 3 sin 2 20i3 sin 2 A 



y = ^ (d 

in the (A, Y) plane, where P and sin 2 A are constant. The intersection of Eq. $1} and 
Y = I/S23 = 2 in the (sin 2 2#i 3 , l/s 23 ) plane is a unique point, which corresponds to a 
solution with eightfold degeneracy. The solution is depicted in Fig.Efa). 



3 



B. cos 2023 + 0, Am^/Am^ = 0, AL = 

At present the Superkamiokande atmospheric neutrino data gives the allowed region 
0.90 < sin 29 23 < 1.0 at 90%CL [jjj, and sin 2 2# 23 can be in general different from 1.0. If 
sin 2 26*23, which is more accurately determined from the oscillation probability P(y M — > v^) 
in the future long baseline experiments, deviates from 1, then we have two solutions for 
Y = 1/4,: 

2 

1 - y/l - sin 2 20 23 
2 

1 + Vl-sin 2 2^ 3 ' 

In this case there are two solutions, the one given by Eq. and Y = Y + and another 
given by Eq. ((TJ and Y = YL. These are two solutions with fourfold degeneracy. The two 
solutions in the (sin 2 26 l 1 3, l/s| 3 ) plane are shown in Fig.^b). From this we see that even 
if we know precisely the values of Piy^ — > v e ), Pip^ — > u e ) and Piy^ — > there are two 
sets of solutions, and this is the ambiguity due to the 6* 2 3 <-> 7r — 9 2 3 degeneracy. 




t_>. cos 26*23 + 0, AmlJAmj-L / 0, AL = 

If we turn on the effect of non-zero Am^ in addition to non-zero 00326*23, then the 
oscillation probabilities are 2 



x 2 sin 2 A + 2xyA sin A cos (5 ± A) + y 2 A 2 



which are correct to the second order in the small parameters lAm^/Am 2 ^ and sin 26?i 3 , 
where 



x = s 2 3 sin 28 



13 



Am 2 x 



Am 31 



c 23 sin 26*12. 



(2) 



In this case, the trajectory of Piy^ 



P, P{y, 



P, where P and P are 



constant, in the (X = sin 29 13, Y = l/s| 3 ) plane is given by a quadratic curve: 



16C X(Y - 1)A sin 2 A 



{P -PfY 2 



sin 2 A 



[(P + P - 2C )(Y - 1) + P + P - 2X sin 2 A] 2 , (3) 



cos 2 A 



where 



Cn 



2 This is obtained by taking the limit A = V2G F N e -> in Eq. (16) in Ref. 0- 
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Eq. © becomes a hyperbola for most of the range of A, but it becomes an ellipse for some 
region A ~ tt. 

When sin 2 2^23 = 1, there are two solutions for the intersection of Y = 2 and Eq. Q. 
This indicates that even if we know the precise values of Piv^ — > u e ), Pivp, — > v e ) and 
P( u ii —> there are two sets of solutions for (#13, #23, with fourfold degeneracy when 
sin 2 2^23 = 1, as is depicted in Fig.|Hfa). This is the ambiguity due to the intrinsic (9\ 3 , 8) 
degeneracy When sin 2 26*23 7^ 1, there are four sets of solutions with twofold degeneracy as 
is depicted in Fig.^b). 



D. cos 2023 + 0, Am^/Aml], / 0, AL ^ 

Furthermore, if we turn on the matter effect AL, then the oscillation probabilities are 
given by 0, [H| 



x 2 f 2 + 2xyfg cos (5 + A) + y 2 g 2 



2„2 



Pip* -> *e) = x 2 / 2 + 2^/(7 cos (5 - A) + y 2 (7 



(4) 



for the normal hierarchy while 



P{v„ 



x 2 f 2 - 2xy]g cos (5 - A) + yV 
x 2 / 2 - 2xyfg cos (5 + A) + y 2 g 2 . 



for the inverted hierarchy where x and y are given by Eq. (J2J), and 



sin (A =f AL/2) 
(1 =f AL/2A) 

sin (AL/2) 
AL/2 A 



(5) 

(6) 
(7) 



Eqs. (HJ) and © are correct up to the second order in lAm^/Am^l and sin 2# 13 , and all 



orders in AL. The trajectory of Pwn 



P, where P and P are 



constant, in the (X = sin 2 29 13 , Y = l/s 2 3 ) plane is again a quadratic curve for either of the 
mass hierarchies: 



IQCXiY 



cos 2 A 
1 



sin 2 A 



P-C P-C 
P-C P-C 



p p n2 



/ 



(Y-l)-(f + f)X + j + j 

(y - 1) - (/ - f)x + 



-1 2 



for the normal hierarchy, and 



16CX(y 



cos 2 A 
1 

sin 2 A 



P-C P-C 



f f 
P-C P-C 



(Y-l)-(f + f)X + j + j 



f 



f 



P P n2 



(Y-l)-(f-f)X + j-j 



(9) 
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for the inverted hierarchy, where 



C 



(Am 2 21 s 


2 

) 


"sin(AL/2)~ 


\Am 2 31/ 




AL/2A 



2 



sin 2 2# 12 . (10) 



Again these quadratic curves become hyperbolas for most of the region of A, but they 
become ellipses for some A ~ n. 

If sin 2 2#23 = 1, then there are four solutions with twofold degeneracy, as is shown in 
Fig.|3^a). If we know for some reason (e.g., from reactor experiments) which solution is 
selected for each mass hierarchy, there are only two solutions. This is the ambiguity due 
to the sgn(Am 31 ) degeneracy. If sin 2 2# 2 3 7^ 1 and if we do not know which solution is 
favored with respect to the intrinsic degeneracy for each hierarchy, and if we do not know 
sgn(Am| 1 ), then there are eight solutions without any degeneracy, as is depicted in Fig.EJb). 
The advantage of our plot is that all the eight solutions for (#13,023) give different points, 
and all the lines in the (sin 2 26 , 13 , l/sl 3 ) plane are described by (at most) quadratic curves 
so that their behaviors are easy to see. 

E. Oscillation maximum 

Finally, let us consider the case where experiments are done at the oscillation maxi- 
mum, i.e., when the neutrino energyi? satisfies A = \Am 3l \L/4E = tt/2. In this case, the 
probabilities become 

Piy^Pe) = x 2 f 2 -2xyfgsm5 + y 2 g 2 , (11) 
Pip^Pe) = x 2 f 2 + 2xyfgsm6 + y 2 g 2 } (12) 

for the normal hierarchy, and 

PW - v e ) = x 2 ] 2 - 2xy]g sin 6 + y 2 g 2 , (13) 
Pfa - Ve) = x 2 f 2 + 2xyfg sin 5 + y 2 g 2 , (14) 

for the inverted hierarchy, where x and y are given by Eq. Q, and f,f,gm Eqs. ©, (0) 
become 

f\_ ^cos(AL/2) ^ _ sin(AL/2) 



/ J lTAL/n ' y AL/n 

for A = 7r/2. The trajectory of Piy^ — > u e ) = P, P{p^ — > v e ) = P in the (X = sin 2 26 , 13 , 
Y = l/s 2 3 ) plane becomes a straight line and is given by 

Y = _/ + / ( y C _) (ik) 

p/f + p/f - c(i/f + i/f) V ff) y ' 

for the normal hierarchy, and 

Y= jj l£ 7tc(x-£) (16) 



p/f + p/f-c(i/f + i/f) V // 

for the inverted hierarchy, where C is given in Eq. ([10)1 . The straight lines (fTB^) and 
(fT6*|) are extremely close to each other in relatively short long baseline experiments such 
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as JPARC, where the matter effect is small. As is shown in Appendix [Bj f|15|) and (|16|) 
have the minimum values in Y = l/s| 3 which is larger than the naive value 1 for either of 
the mass hierarchies. Since Eqs. (115)1 and (|16|1 are linear in X, there is only one solution 
between them and F=const. Thus the ambiguity due to the intrinsic degeneracy is solved 
by performing experiments at the oscillation maximum, although it is then transformed into 
another ambiguity due to the 5 «-> n — 5 degeneracy. 

If sin 2 2#23 — 1, then all the four solutions are basically close to each other in the (sin 2 2# 13 , 
I/S23) plane, and the ambiguity due to degeneracies are not serious as far as #13 and #23 are 
concerned (See Fig.^Ja)). On the other hand, if sin 2 2 #23 deviates fairly from 1, then the 
solutions are separated into two groups, those for #23 > tt/4 and those for #23 < vr/4 in the 
(sin 2 2#i3, I/S23) plane, as is shown in Fig.|Jb). In this case resolution of the #23 *-> n/2 — 623 
ambiguity is necessary to determine #13, #23 and S. 



F. Fake effects on CP violation due to degeneracies 

1. sin 2 26»23~l 

If the JPARC experiment finds out from the measurement of the disappearance proba- 
bility Piy^ — > v^) = P that sin 2 26* 2 3 — 1.0 with a good approximation, then we would not 
have to worry very much about parameter degeneracy as far as #13 and #23 are concerned, 
since the values of #13 and #23 for all the different solutions are close to each other. 

On the other hand, when it comes to the value of the CP phase phase 5, we have to be 
careful. From Ref. [9| the true value 5 and the fake value 5' for the CP phase satisfy the 
following: 

, • w • s f + P-ff -r 2 P + P f-f 

x smo = xsmo = - A = , 17 

// sin A // 2yg ' 1 } 

where x, y are given in Eq. ((2)), /, /, g are given in Eqs. (JSJ) and (jZJ), and x' is defined by 

x a = x\P + P ~ ff) ~ 2y_g(f - f)x sin 5 sin A 

Eq. (|17p indicates that even if sin S = we have nonvanishing fake CP violating effect 



sin^ = -x /2+ / 2 f L J H (18) 

// 2^sinA\// 2 + / 2 -//' 1 ^ 

if we fail to identify the correct sign of Am^. In the case of the JPARC experiment, Eq. 
(fTHj) implies 



sin 5' ~ — 2.2 sin 



13) 



which is not negligible unless sin 2 2# 13 <C 10 2 . Therefore we have to know the sign of Am^ 
to determine the CP phase to good precision. 
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2. sin 2 2023 < 1 



As was explained in Sect. Ill El if sin 2 2023 deviates fairly from 1, then we have to resolve 
the ambiguity due to the 2 3 degeneracy to determine the values of 6 Vi and 2 3 • As for the 
value of the CP phase 5, we can estimate how serious the effect of the 023 ambiguity on the 
value of 5 could be. If the true value 5 is zero, then the CP phase 5' for the fake solution 
can be estimated as [9| 



sin 26^3 sin 5' 



A77-4 



Am 31 



sin 20i2 cot 20 23 



sin A 



where 



sin 2 20^ = sin 2 20! 3 tan 2 



13 



'23 



a?t4 

Am 2 x 



g 2 sin 2 20 



12 



1 - tan 2 23 ) 



and /, /, g are defined in Eqs. © and (JJJ). In the case of JPARC, we have 



I sin 5' I 



1 I cot 20 



23 



1 



< 



1 



1 



200 t 23 sin20 13 500 J^29 



(19) 



13 



where we have used the bound 0.90 < sin 2 2023 < 1.0 from the atmospheric neutrino data 
in the second inequality, so that we see that the ambiguity due to the 023 does not cause a 
serious problem on determination of S for sin 2 20 13 > 10~ 2 . It should be stressed, however, 
that the effect on CP violation due to the sgn(Am| 1 ) ambiguity is also serious in this case. 



III. RESOLUTION OF AMBIGUITIES BY THE THIRD MEASUREMENT AF- 
TER JPARC 

In this section, assuming that the JPARC experiment, which is expected to be the first 
superbeam experiment, measures P(u^ — > v e ) and P(pn — > v e ) a ^ the oscillation maximum 
A = Am^L/AE = tt/2, we will discuss how the third measurement after JPARC can 
resolve the ambiguities by using the plot in the (sin 2 20 13 , l/s 2 3 ) plane. Resolution of the 
023 ambiguity has been discussed using the disappearance measurement of P{p e — > u e ) at 
reactors jaEHIliJ 12, 3 and the silver channel P{y e — > v T ) at neutrino factories [Hj, but 
it has not been discussed much using the channel —>■ v e 3 . Here we take the following 
reference values for the oscillation parameters: 

sin 2 20i2 = 0.8, sin 2 20 i3 = 0.05, sin 2 20 23 = 0.96, 

Am 2 x = 7 x 10" 5 eV 2 , Am 2 x = 2.5 x 10~ 3 eV 2 > 0, 5 = vr/4. (20) 



3 There have been a lot of works |2fll | on how to resolve parameter degeneracies, but they discussed mainly 
the intrinsic and sgn(Am| 1 ) degeneracies, and the present scenario, in which the third experiment follows 
the JPARC results on P(f M — » v e ) plus P(y^ — > v e ) which are measured at the oscillation maximum, has 
not been considered. 
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A. -> v e 



Let us discuss the case in which another long baseline experiment measures P{p 11 — > v e )- 
From the measurements of Piv^ — > f e ) and Pip^ — > u e ) by JPARC at the oscillation 
maximum we can deduce the value of 5, up to the eightfold ambiguity (5 <-> n — 5, 9 2 3 *-> 
7r/2 — 623, Am^j <-> -Am^) 4 As is depicted in Fig.|BJ depending on whether s 2 ^ — 1/2 
is positive or negative, we assign the subscript ±, and depending on whether our ansatz 
for sgn(Amg 1 ) is correct or wrong, we assign the subscript c or w. Thus the eight possible 
values of 5 are given by 

S +w , 5 +c , <5_ w , 5-c, 71 - S +w , 71 - S +c , 7T - <L W , 7T - 5- c . (21) 

Now suppose that the third measurement gives the value P for the oscillation probability 
P( u ii ~^ v e)- Then there are in general eight lines in the [X = sin 2 29 13, Y = l/s| 3 ) plane 
given by 

fX = [P — C + 2C cos 2 (5 + A)] (Y — 1) + P 



-2 cos(5 + A) y/C(Y -l)yJ[P-C sin 2 (5 + A)] (Y - 1) + P (22) 
for the normal hierarchy, and 

] 2 X = [P-C + 2Ccos 2 {5 - A)] (y - 1) + P 



-2 cos((5 - A) y/C(Y -l)yJ[P-C sin 2 (5 - A)] (Y - 1) + P (23) 

for the inverted hierarchy, where C is defined in Eq. ([10|h A = |Am 2 1 |L/4i? is defined for 
the third measurement, and 5 takes one of the eight values given in Eq. (j21j) . The derivation 
of and is given in Appendix |X| It turns out that the solutions and are 
hyperbola if cos 2 (5 ± A) > (C — P)/P, where + and — refer to the normal and inverted 
hierarchy, and ellipses if cos 2 (5 ± A) < (C — P)/P. In practice, however, the difference 
between hyperbola and ellipses is not so important for the present discussions, because we 
are only interested in the behaviors of these curves in the region 1.52 < Y = l/s| 3 < 2.92 
which comes from the 90%CL allowed region of the Superkamiokande atmospheric neutrino 
data for sin 2 26*23 . 

Here let us look at three typical cases: L=295km, L=730km, L=3000km, each of which 
corresponds to JPARC, off-axis NuMI [2^, and a neutrino factory 2^1 5 



4 I thank Hiroaki Sugiyama for pointing this out to me. 

5 For L=3000km the density of the matter may not be treated as constant, and the probability formulae 
and (JSJ may no longer be valid. It turns out, however, that the approximation of the formulae becomes 
good if we replace AL by AL — > J Q L A{x)dx everywhere in the formula. In the following discussions, the 
replacement AL — > J Q A(x)dx is always understood in the case of the baseline L=3000km. It should be 
mentioned that the neutrino energy spectrum at neutrino factories is continuous and it is assumed here 
that we take one particular energy bin whose energy range can be made relatively small. It should be 
also noted that neutrino factories actually measure the probabilities P{v e — > v^) or P(y e — > v^), instead 
of P(Vfj, — * 24) or P(Vfj, — » v e ). Here we discuss for simplicity the trajectory of Piy^ — > v e ) whose feature 
is the same as that of P(v e — » v^). 



9 



Figs. 1718191 show the trajectories of P{p IJj — > v e ) obtained in the third measurement to- 
gether with the constraint of Piy^ — > is e ), Pipy, — > v e ) and Piy^ — > v^) by JPARC, for 
L=295km, L=730km, L=3000km, respectively, where A = \/S.m\ 1 \L / AE takes the values 
A = jn/8 (j = 1, • • • ,7). The purple (light blue) blob stands for the true (fake) solution 
given by the JPARC results on -P(z^ — > u e ), P{p^ — > v e ) an d P{ u n ~> u ^)- For the correct 
(wrong) guess on the mass hierarchy, there are in general four red (blue) curves because 
the CP phase 5, which is deduced from the JPARC results on P{u^ — > u e ), P{p^ — > u e ) 
and Piy^ — > Vp), is fourfold: (S +c , 5- c , tt — 5 +c , tt — <5_ c ) for the correct assumption on the 
hierarchy and (5+ w , <5_ w , tt — 5 +w , 7r — 5_ w ) for the wrong assumption. In most cases the four 
(red or blue) curves are separated into two pairs of curves. As we will see later, the large 
split is due to the 5 <-» tt — 5 ambiguity, while the small split is due to the 6> 2 3 <-> tt/2 — 6> 2 3 
ambiguity. The reason that the latter splitting is small is because the difference of the values 
in the CP phases is small, as is seen from Eq. (119)) . In some of the figures in Figs. 171819*1 the 
number of the red or blue curves is less than four because not all values of S give consistent 
solutions for a set of the oscillation parameters. 

Let us discuss each ambiguity one by one. 

1. 5 <— > 7r — 5 ambiguity 

As was mentioned above, the large splitting of four (red or blue) lines into two pair of lines 
is due to the 5 <->• tt — 5 ambiguity. From Eqs. (|22|) and (|23j) we see that the only difference 
of the solutions with 5 and with tt — 5 appears in cos(S ± A) or sin(<5 ± A). If A = tt/2 (i.e., 
the oscillation maximum), we have cos(<5 + A) = — sin 5, cos(7r — 5 + A) = — sin 5, so that the 
values of X with 5 and with tt — 5 are the same, i.e., at oscillation maximum there is exact 
5 «-> tt — 6 degeneracy. On the other hand, if A 7^ tt/2, we have cos(<5 + A) 7^ cos(7r — 5 + A), 
and the values of X with S and with tt — 5 are different. Thus, to resolve the 5 «-> tt — 5 
ambiguity it is advantageous to perform an experiment at A which is farther away from 
tt/2. Deviation of A from tt/2 implies either high energy or low energy. In general the 
number of events increases for high energy because both the cross section and the neutrino 
flux increase, so the high energy option is preferred to resolve the 5 «-> tt — 5 ambiguity 6 . 

2. Amjj «-> -Ara^ ambiguity 

As one can easily imagine, the sgn(Amg 1 ) ambiguity is resolved better with longer base- 
lines, since the dimensionless quantity AL = \f2GpN e L ~ (L/1900km)(p/2.7g-cm -3 ) be- 
comes of order one for L ->1000km. On the other hand, from Figs.|Hl and OH we observe 
that the split of the curves with the different mass hierarchies (the red vs blue curves) is 
larger for lower energy. Naively this appears to be counterintuitive, because at low energy 
the matter effect is expected to be less important (|Am§ 1 |L/4E ^> AL). However, this is 
not the case because we are dealing with the value of sin 2 2#i3 which is obtained for a given 
value of P(Vfj, — > v e ). To see this, let us consider for simplicity the the value of X = sin 2 2#i3 
at Y = l/s 2 3 — 1) i- e -) the X-intercept of the quadratic curves at Y — 1. (sin 2 2# 13 ) n 
((sin 2 26 , 13 ) i ) at Y — 1 for the normal (inverted) hierarchy is given by x 2 by putting y = 



6 Resolution of 6 <-> n — S ambiguity at neutrino factories was discussed in [l 
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in Eq. (gj) (Eq. ©): 



(sin 2 2# 13 ) n = j 2 



for S23 = 1 



(sin 2 20 13 )i = j r 
The ratio of these two quantities is given for small AL by 



^sin 2 2# 13 ) n / 2 _ sin 2 (A-AL/2) /1 + AL/2AY 
(sin 2 2# 13 ), ~ T 2 ~ sin 2 (A + AL/2) \l - AL/2A ) 



, 1 1 

1 + 2AL 



A tan A 

so that the larger A is (the smaller the neutrino energy is), the larger this ratio becomes, 
as long as A does not exceed 7r. This phenomenon suggests that it is potentially possible to 
enhance the matter effect by performing an experiment at low energy (A > n/2) even with 
L=730km, and it may enable us to determine the sign of Am 31 at the off-axis NuMI exper- 
iment. While the neutrino flux decreases for low energy at the off-axis NuMI experiment, 
the cross section at E ~lGeV is not particularly small compared to higher energy, so the 
low energy possibility at the off-axis NuMI experiment deserves serious study 



3. #23 <-> — #23 ambiguity 

Figs. OHM which are plotted for 5 = 7r/4, suggest that there is a tendency in which the 
slope of the red curve which goes through the true point (the purple blob) is almost the 
same for high energy as that of the straight green line obtained by JPARC, while for the 
low energy the slope of the red curve is smaller than that of the JPARC green line. Here 
we will discuss the X-intercept at Y = 1 instead of calculating the slope itself, because it is 
easier to consider the X-intercept and because the difference in the X-intercepts inevitably 
implies the different slopes for the two lines, as almost all the curves are approximately 
straight lines. In the case of JPARC, the matter effect is small (AL ~ 0.08) so that we can 
put / ~ / ~ 1. From Eq. (|15J1 we have the X-intercept at Y — 1 

XjFARC = p Jl±^l s f+f = 1 ., (24) 

where the term g 2 y 2 has been ignored for simplicity. On the other hand, for the third 
measurement, from Eq. (|22|) we have 

X 3rd = ~x 2 + 2^xycos(5 + A), (25) 

where the term g 2 y 2 has been ignored again for simplicity. Eq. (J25)) indicates that it is the 
second term in Eq. (J2*5j) that deviates the intercept X 3r< j of the red line from the intercept 
Xjparc of the JPARC green line. In order for the difference between Xjparc and X 3r d to 
be large, / has to be small and | cos(<5 + A)| has to be large. When AL is small, in order for 
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/ to be small, WAm^L/AE — ir\ has to be small. This is one of the conditions to resolve 
the 823 ambiguity. Here we are using the reference value 5 = tt/4, so the deviation becomes 
maximal if \8 + A\ = |7r/4 + A| ~7r. In real experiments, however, nobody knows the value 
of the true 6 in advance, so it is difficult to design a long baseline experiment to resolve the 
#23 ^ tt/2 — #23 ambiguity. If 5 turns out to satisfy | cos(5 + A)| ~ 1 in the result of the 
third experiment, then we may be able to resolve the 623 ambiguity as a byproduct. 



B. 



It turns out that the situation does not change very much even if we use the 



channel in the third experiment. Typical curves are given for z/„ 



in Fig.^J which are 



similar to those in Fig. l7l8liJl Thus the conclusions drawn on resolution of the ambiguities 
hold qualitatively in the case of z/„ — > v e channel. 



The experiment with the channel v e — > v T requires intense v e beams and it is expected 
that such measurements can be done at neutrino factories or at beta beam experiments I2F 
The oscillation probability P[y e — > v T ) is given by 



P{u e -> v T ) = x 2 f 2 + 2fgxy cos(5 + A) 



~2 2 

y 9 , 



where 



x 



y 



c 2 3 sin 26*23 



Am 2 21 



S23 sin 29 



12; 



and /, g are given in Eqs. (JHJ) and Q. The solution for P{y e 
constant, is given by 



u T ) = Q, where Q is 



X 



Q 
f 2 



1 + 



2cos 2 (^ + A) 
1-C/Q 

2 cos(5 + A) , 

' 1 + 



1-C/Q 
Y - 1 



+ 1 



1-C/Q 



cos 2 (5 + A) 
1-C/Q 



C/Q 



Y - 1 



+ 1 



(26) 



where X = sin 2 26*13, Y = l/ s 23 as before and C is given in Eq. (fTUj) . Eqs. is plotted 
in Fig. ^2 i n the case of L=2810km. From Fig.E] we see that the curve P{y e — > v T ) = 
Q intersects with the JPARC green line almost perpendicularly, and it is experimentally 
advantageous. Namely, in real experiments all the measured quantities have errors and the 
curves become thick. In this case the allowed region is small area around the true solution 
in the (sin 2 26*13, l/ s 23) plane and one expects that the fake solution with respect to the #23 
ambiguity can be excluded. This is in contrast to the case of the —>■ v e and — > v e 
channels, in which the slope of the red curves is almost the same as that of the JPARC 
green line and the allowed region can easily contain both the true and fake solutions, so that 
it becomes difficult to distinguish the true point from the fake one. 
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As in the case of the — > v e channel, the 5 <-> tt — 5 ambiguity is expected to be resolved 
more likely for the larger value of | A — n/2\, and the sgn(Am| 1 ) ambiguity is resolved easily 
for larger baseline L (e.g., L ~3000km). 

Thus the measurement of the v e — > v T channel is a promising possibility as a potentially 
powerful candidate to resolve parameter degeneracies in the future. 

IV. DISCUSSION AND CONCLUSION 

In this paper we have shown that the eightfold parameter degeneracy in neutrino oscilla- 
tions can be easily seen by plotting the trajectory of constant probabilities in the (sin 2 2# 13 , 
1/533) plane. Using this plot, we have seen that the third measurement after the JPARC 
results on P{u^ — > u e ) and P(u^ — > v e ) may resolve the sgn(Amg 1 ) ambiguity at L >1000km, 
the 5 <-> 71 — 5 ambiguity off the oscillation maximum (|A — 7r/2|~0(l)), and the 623 am- 
biguity if ||Am| 1 |L/4i? — 7r I is small and S turns out to satisfy |cos(<5 + A)| ~ 1. In general 
all these constraints on A = |Am| 1 |L/4E may be satisfied by taking A = 71. The condi- 
tion A = 7T, however, actually corresponds to the oscillation minimum, and the number of 
events is expected to be small for a number of reasons: (1) The probability itself is small 
at the oscillation minimum; (2) A = n implies low energy and the neutrino flux decreases 
at low energy; (3) The cross section is in general smaller at low energy than that at high 
energy. Therefore, to gain statistics, it is presumably wise to perform an experiment at 
tt/2 < A < 71 after JPARC. The off-axis NuMI experiment with tt/2 < A < tt (E ~lGeV) 
may have advantage to resolve these ambiguities. 

As is seen in Figs. |B1 and the experiments at the oscillation maximum does not appear 
to be useful after JPARC except for the sgn(Am| 1 ) ambiguity. In order to achieve other 
goals such as resolution of the 5 <-> n — 5 ambiguity and the #23 ambiguity, it is wise to stay 
away from A = 71/ 2 in experiments after JPARC. 

Although only the oscillation probabilities were discussed without taking the statistical 
and systematic errors into account in this paper, we hope that the present work gives some 
insight on how the ambiguities may be resolved in the future long baseline experiments. 

APPENDIX A: EXPRESSION FOR P(z^ u e ) = P 

First of all, let us derive Eqs. (|22|) and (|23|). For the normal hierarchy, the probability 
P{vn -» f e ) = P is given by 

P = P( Vfl -> v e ) = x 2 f + 2xyfg cos (5 + A) + y 2 g 2 

= f 2 y + 2f S\f^-y) cos {6 + A) + c " f) ' (A1) 

where X = sin 2 29\3, Y = l/s| 3 as in the text, / is defined in Eq. (JBJ), and C is given in Eq. 
(JTUJ). Eq. ()A1|) is rewritten as 

(P - C){Y - 1) + P - f 2 X = 2^f 2 ~X^C(Y - l)cos (5 + A) . (A2) 
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Taking the square of the both hand sides of Eq. (|A2[) . we get 7 

{fX) 2 - ifX { [P - C + 2Ccos 2 {5 + A)] (Y - 1) + P} + [(P - C)(Y - 1) + P} 2 = 0. 

Solving this quadratic equation, we obtain 

f 2 X = [P-C + 2Ccos 2 {5 + A)] (Y - 1) + P 

±^/{[P - C + 2CW (5 + A)] (Y — 1) + P} — [(P - C)(Y - 1) + P] 2 
= [P - C + 2Ccos 2 (5 + A)] (Y - 1) + P 

±2 cos (5 + A) y/C(Y -l)y/[P-C + C cos 2 {5 + A)] (Y - 1) + P. (A3) 

If cos (5 + A) > then from Eq. {£2J we see that (P + C)(Y - 1) + P - / 2 X has to be 
positive. On the other hand, Eq. (jA3|) gives 



(P + C)(Y -1)+P- fX 
T2 cos (5 + A) y/C(Y - 1) 

'[P - C + Ccos 2 (5 + A)] (Y - 1) + P - cos {5 + A) y/C{Y-l)} . (A4) 



From Eq. (|A4|) we conclude that we have to take the minus sign in Eq. (|A3|) for the right 
hand side of Eq. (jA4j) to be positive. Hence from Pivp, — ► u e ) = P we get 

fX = [P-C + 2Ccos 2 (5 + A)] (Y-l) + P 



-2 cos(5 + A) y/C(Y - 1)J[P-C sin 2 (5 + A)] (Y - 1) + P, (A5) 



and from P{y e — > v^) = P we have 



fX = [P-C + 2Ccos 2 (5 - A)] (Y - 1) + P 



-2 cos(<5 - A) v/C(Y - l)yf [P - C sin 2 (5 - A)] (Y - 1) + P. (A6) 

When P — C sin 2 (5 + A) > 0, Eq. ()A5|) is a hyperbola, and the physical region for Y — 1 is 
Y — 1 > 0. On the other hand, when P — C sin 2 (5 + A) < 0, Eq. (|A5|) becomes an ellipse 
and the physical region for Y — lisO<Y — 1 < P/[C sin 2 (5 + A) — P]. 
Similarly, we obtain for the inverted hierarchy: 

fX = [P-C + 2Ccos 2 (5 - A)] (Y - 1) + P 



+2 cos(<5 - A) y/C(Y-l)yf [P-C sin 2 (5 - A)] (Y - 1) + P, (A7) 
fX = [P-C + 2Ccos 2 (5 + A)] (Y - 1) + P 



+2 cos((5 + A) y/C(Y-l)J [P-C sin 2 ((5 + A)] (Y - 1) + P. (A8) 



7 Here we consider for simplicity the case where all the arguments of the square root are positive. After we 
obtain the final result, we see that the final formula makes sense as long as the whole product of all the 
arguments is positive. 
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APPENDIX B: TRAJECTORIES AT THE OSCILLATION MAXIMUM 



Throughout this appendix we will assume A = n/2 and we will assume Am^ > for 
most part of this appendix. From Eq. (jA5|) . the condition P[y^ — > v e ) = P for the neutrino 
mode alone gives 

f 2 X = [P-C + 2Csin 2 6]{Y-1) + P 

+2 sin 5 y/C(Y -l)y/[P-C cos 2 5} (Y - 1) + P, (Bl) 
while the condition Piv^ — > z/ e ) = P for the anti-neutrino mode alone gives 
] 2 X = [P-C + 2Csin 2 5] (Y - 1) + P 



-2 sin 5 ^C(y ^T) J [P - C cos 2 5] (Y - 1) + P. (B2) 

When 5 ranges from — 7r/2 to 7r/2, Eq. (jBljl sweeps out the inside of a hyperbola, as 
is depicted by the red curves in Fig. H2T a). while ()B2|) sweeps out the inside of another 
hyperbola for the anti- neutrino mode (cf. the blue curves in Fig. lT2T a)). Notice that the 
left (right) edge of the hyperbola (jBl|) for the neutrino mode corresponds to 5 = —tt/2 
(5 = +vr/2) whereas the left (right) edge of the other hyperbola (|B2|) for the anti- neutrino 
mode corresponds to 5 = +tc/2 (5 = —tt/2). Since the straight line ()15|) is the intersection 
of the two regions (the yellow and light blue regions in Fig. lT2T b)). the lowest point in the 
straight line is obtained by putting S = +tt/2 (5 = -tt/2) if P/f 2 < P/f 2 (if P/f 2 > P/f 2 ), 
respectively, depending on whether the region for the anti-neutrino mode is to the right of 
that for the neutrino mode. Therefore, if P/f 2 < P/f 2 , then putting 5 = +tt/2 in Eqs. 
and (fT2")l and assuming xf > yg, which should hold if sin 2 2#i3 is not so small, we get 




which lead to the minimum value of Y 



Y, 



(n) 



(fVP-fVPf 
c(f + f) 2 



for the normal hierarchy. On the other hand, for the inverted hierarchy, the corresponding 
values of 5 for the edges for the two modes are the same as those for the normal hierarchy 
(5 = ±tt/2). Hence, if P/f 2 < P/f 2 , then putting 6 = +vr/2 in Eqs. (JH and (HJ) and 
assuming xf > yg, we obtain 

VP = xf-yg 
\fp = x f + yg 
which leads to the minimum value of Y 



Y 



(i) 



(fVP - fVP) 
C(f + f) 2 



for Amjj < 0. 
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FIG. 1: Two mass patterns, (a), (b) correspond to the normal (Am^j > 0), inverted (A 
hierarchy, respectively. 




FIG. 2: The solutions, which are marked by black blobs, for given P{ya — ► v e ), P(pu & e ) 
and P{Vfj, — > Vn) in the case of Am^/Am^ = 0, AL = 0. (a) For cos 2#23 = 0, the intersection of 
Y = l/s| 3 = 2 and the trajectory of Piy^ — > u e ) = Piy^ — > z^) = const, is one point with eightfold 
degeneracy, (b) For cos 2#23 7^ 0, the intersections are two solutions with fourfold degeneracy. 
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FIG. 3: The solutions, which are marked by black blobs, for given Piy^ — * u e ), Pip^ — * v e ) and 
P( u n —> u n) in the case of Am^/Am^ ^ 0,AL = 0. (a) For cos2#23 = 0, the intersection of 
Y = 1/^23 = 2 and the trajectory of P — ► u e ) = const., -P(f^ — ► Vn) = const, are two points 
with fourfold degeneracy, (b) For cos 26*23 7^ 0, the intersections are four solutions with twofold 
degeneracy. 
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FIG. 4: The solutions, which are marked by black blobs, for given Piy^ — * u e ), Pip^ — * v e ) and 
P( u n —> u n) in the case of Am^/Am^ ^ 0,AL ^ 0. (a) For cos2#23 = 0, the intersection of 
Y = 1/^23 = 2 and the trajectory of P{yu — > v e ) = const., P {yn — ► f„) = const, four points with 
twofold degeneracy, (b) For cos 26*23 0> the intersections are eight solutions without degeneracy. 
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FIG. 5: The #23 ambiguity which could arise after the JPARC measurements of P{v^ — > v e ), 
P(pn — > v e ) and P{vn — * fa) at the oscillation maximum, (a) If sin 2 2^23 — 1.0 then the values 
of 0i3 and #23 are close to each other for all the solutions, and the ambiguity is not serious as far 
as #13 and #23 are concerned, (b) If sin 2 26*23 < 1 then the #23 ambiguity has to be resolved to 
determine #13 and #23 to good precision. 
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FIG. 6: Four possible values for the CP phase 6 at the oscillation maximum. The red (blue) line 
stands for the normal (inverted) hierarchy. Since we are assuming the normal hierarchy here, the 
red (blue) line corresponds to the correct (wrong) assumption on the mass hierarchy, db sign stands 
for the choice of S23 = (1 i yl — sin 2 2^23) /2 in the #23 ambiguity, and c (w) stands for the correct 
(wrong) assumption on the mass hierarchy. 
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L = 295 km jparc v+ v 

true solution • correct 

fake solution • wrong 

FIG. 7: The trajectories of P(vu v e) = const, of the third experiment at L=295km with 
A = lAmj^L/AE = (j/8)ir (0 < j < 7,j + 4) after JPARC. The true values are those in Eq. flUJ). 
The green line is the JPARC result obtained by P{v^ — * v e ) and Piy^ — * v e ) at the oscillation 
maximum. The red (blue) lines are the trajectories of Piy^ — > v e ) given by the third experiment 
assuming the normal (inverted) hierarchy, where S takes four values for each mass hierarchy. 
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FIG. 8: The trajectories of P (i/„ — * v e ) = const, of the third experiment at L=730km with 
A = \Am 2 31 \L/4E = (j'/8)vr (0 < j < 7) after JPARC. The true values are those in Eq. (EOJ). 
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FIG. 9: The trajectories of -P(z^ — ► v e ) = const, of the third experiment at L=3000km with 
A = lAml^L/AE = (j/8)n (0 < j < 7) after JPARC. The true values are those in Eq. (JUJ). For 
A > (3/8)7r, the blue curves (with the wrong assumption for the mass hierarchy) are not in the 
figure because they are far to the right. 
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FIG. 10: The trajectories oiPiy^ — > u e ) = const, of the third experiment with A = \Am^\L/4:E = 
tt/8 after JPARC. The behaviors are almost similar to those for P(va — * v e ) = const. The true 
values are those in Eq. (|2U|). 
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FIG. 11: The trajectories of P{v e — ► v T ) = const, of the third experiment at L=2810km with 
A = \£rnil x \L/AE = (j/8)vr (j = 1,2,3) after JPARC. The true values are those in Eq. The 
curves intersect with the JPARC line perpendicularly, so this channel is advantageous to resolve 
the ambiguities from experimental point of view. 
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FIG. 12: The region of constant probabilities at the oscillation maximum, (a) Each red (blue) line 
stands for P(vn — ¥ v e) = const. [P{pn v e ) = const.) for a specific value of 5. The red line on the 
right (left) edge corresponds to S = +n/2 (<5 = — vr/2), while the blue line on the edge right (left) 
corresponds to 6 = —ir/2 (5 = +ir/2). (b) When 5 varies from to 2n, the line — > f e ) = 
const, sweeps out the yellow region, whereas the line — > D e ) = const, sweeps out the light blue 
region. The black straight line, which is given by P{v^ — > v e ) = const, and Pip^ — > v e ) = const., 
lies in the overlapping green region. 



